Magnetic Component of Quark-Gluon Plasma is also a Liquid! 
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The magnetic scenario recently suggested in [1] emphasizes the role of monopoles in strongly coupled 
quark-gluon plasma (sQGP) near/above the deconfinement temperature, and specifically predicts 
that they help reduce its viscosity by the "magnetic bottle" effect. Arguments for "magnetic liquid" 
in l-2Tc based on lattice results of monopole density were provided in [2]. Here we present results for 
monopole-(anti)monopole correlation functions from the same classical molecular dynamics simula- 
tions, which are found to be in good agreement with lattice results in [3] . We show that the magnetic 
Coulomb coupling runs in the direction opposite to the electric one, as expected, and it is roughly 
inverse of the asymptotic freedom formula for the electric one. However, as T decreases to Tc, the 
magnetic coupling never gets too weak, with the plasma parameter always large enough (F > 1). 
This nicely agrees with empirical evidences from RHIC experiments, implying that magnetic objects 
cannot have large mean free path and should also form a good liquid with low viscosity. 
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I. RUNNING COUPLINGS 

Discussions of electric-magnetic duality appear in the- 
oretical literature regularly since Maxwell's time. Espe- 
cially important for what follows are: (i) the celebrated 
Dirac quantization condition [4], (ii) 't Hooft-Polyakov 
monopole solutions [7], (iii) "dual superconductor" idea 
of confinement by 't Hooft and Mandelstam and (iv) 
Seiberg-Witten solution of the J\f = 2 SUSY gauge theory 
[6] , identifying properties and dynamical role of magnet- 
ically charged objects in this setting. 

Before we focus on our main subject - monopole mu- 
tual interactions - let us first comment on other theoreti- 
cal issues related with situations in which both "electric" 
and "magnetic" particles are present at the same time. 

For pure gauge fields electric/magnetic duality simply 
means rewriting magnetic field B as gradient and elec- 
tric field as a curl of a "dual" potential: however there 
are nontrivial questions about the sources (and bound- 
ary conditions). The electric objects - quarks - are tra- 
ditionally present in the Lagrangian as Noether charges, 
while monopoles are solitonic solutions carrying topolog- 
ical charges. Can "magnetic" formulation be consistently 
defined, interchanging their roles and putting monopoles 
in the Lagrangian instead? Can even a situation be found 
in which both formulations are similar? These ideas were 
discussed starting from the famous paper [8]. 

Since both description should describe the same theory 
serious issues of consistency appear. At the quantum me- 
chanics level the famous Dirac [4] condition must be held, 
demanding basically that the product of two couplings is 
fixed to an integer (put here to f ) 



as * aM 



1 



(1) 



So while one of them may be small, the other must nec- 
essarily be large. 



At the level of quantum field theory the Dirac condi- 
tion elevates into a requirement that two couplings must 
run into the opposite directions 



Pia)E + /3(a) Af = 



(2) 



where the beta functions are (3{a) = — 2(3{g)/g 

for the electric and magnetic couplings respectively, with 
f3{g) = being the usual beta function. This indeed is 
what happens in Seiberg-Witten solution, in which elec- 
tric coupling is weak at large momenta due to asymptotic 
freedom, and magnetic is weak at small ones due to U(f ) 
"Landau pole" . 

In our previous paper [1] we discussed for sQGP the 
important role of this generic feature. As it is known for 
30 years, QGP at high T can be described perturbatively, 
with e.g. small quark and gluon effective masses M/T ~ 
^/oieiectric << 1 Thc monopolcs in this case are heavy 
composites which play a minor role, although they are 
strongly interacting and form an interesting sub-sector 
in which perturbative analysis is impossible. However 
as T goes down and one approaches the deconfinement 
transition T ^ Tc^ the inverse is expected to happen: 
electrically charged particles - quarks and gluons ~ are 
getting heavier and more strongly coupled. There are 
strong evidences that both phenomena do happen: lattice 
data show that quark (baryon) masses seem to be large in 
sQGP near Tc [9], while RHIC's "perfect liquid" [10,11] 
supports the idea of strong coupling. 

However at this point one may ask what happens with 
monopoles: as T ^ Tc the same logic suggests that they 
must become lighter and more important. At some point 
their masses (and roles) get comparable to that of the 
electric objects, after which the tables are turned and 
their fortune reversed. Electric objects gets strongly cou- 
pled and complicated while monopoles gets lighter, pro- 
liferate and eventually take over the bulk, expelling elec- 
tric fields into the flux tube. As shown in [12], this may 



happen in the plasma phase, before confinement transi- 
tion is reached. The corresponding phase diagram was 
discussed in [1]. Here are two questions, on which will 
we be focused below: (i) Are there evidences that the 
magnetic coupling does run in the opposite direction? 
(ii) How small does the magnetic coupling become at 
T Tc, and is a perturbative description of magnetic 
plasma possible? As the reader will see, we will answer 
"yes" to the first and "no" to the second question. 

(A digression about the most symmetric Af—4: SYM 
theory which is conformal. How do we know that its cou- 
pling does not run? One may calculate the first coefficient 
of the beta function, and will indeed see that negative 
gauge contribution is nicely cancelled by fermions and 
scalars. But there are infinitely many coefficients, and 
one has to check them all! An elegant way to prove the 
case is based on another outstanding feature of the A/'=4 
SYM: this theory is (nearly) self — dual under electric- 
magnetic duality. As we discussed above, the Dirac con- 
dition requires the product of electric and magnetic cou- 
plings to be constant. But, as shown by Osborn [13], 
in Higgs case the multiplet of (lowest) magnetic objects 
of the Af=4: SYM theory include 5 scalars, 4 fermions 
(monopoles plus one gluino zero mode occupied), plus 1 
spin-1 particle (3 polarizations), exactly the same set of 
states as in the original electric multiplet (giuon-gluinoes- 
Higgses). Thus effective magnetic theory has the same 
Lagrangian as the original electric formulation and the 
same beta function. That would conflict with require- 
ment that both couplings run in the opposite direction, 
unless they do not run at all!) 



II. MOLECULAR DYNAMICS WITH 
MONOPOLES 

It has became apparent around 2003 [10,11] that 
Quark-Gluon Plasma discovered at Relativistic Heavy 
Ion Collision (RHIC) is the most "perfect liquid" known. 
Theorists have since been taking the big challenge to ex- 
plain the remarkable properties of sQGP with diverse 
approaches, ranging from those borrowed from classical 
plasmas to AdS/CFT duality: for recent reviews see [14]. 

Molecular Dynamics (MD) simulations are proved to 
be powerful tools widely used for studying the conven- 
tional EM plasma, especially in the strongly coupled 
regime when the analytic approaches are difficult. It pro- 
vides detailed microscopic real time information about 
correlation functions and transport properties. It has re- 
cently been employed to simulate sQGP in [15,1]. The 
key classical Coulomb plasma parameter F is defined by 



(3) 



to the thermal kinetic energy; thus a weakly coupled 
plasma has F << 1. Plasmas with F > 1 are known 
as strongly coupled, and for F ~ l..Fc ^ 100 it is in the 
liquid regime, becoming solid for F > Fc. Less precisely 
defined "glassy" regime is at Gamma = 10.. F'^ in be- 
tween. The electric F parameter in sQGP is believed to 
be Te ~ 3 at T = (1 — 2)Tc [15], so it corresponds to a 
liquid. The value of the magnetic parameter Tm is the 
subject of this Letter. 
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FIG. 1. (color online) Monopole-antimonopole (boxes with 
red curve) and monopole-monopole (triangles with blue 
curve) correlators from MD simulations at F = 2.3 [1]. 

In our previous paper [1] we explored (to our knowl- 
edge, for the first time) a novel plasmas made of a mix- 
ture of both electric and magnetic charges. Using MD 
and standard Kubo formulae we calculated its transport 
properties, i.e. the shear viscosity and diffusion constant. 
We found that an equal 50%-50% mixture has smallest 
viscosity (the shortest mean free path, if one want to use 
such language) because of the "magnetic bottle' effect. 

Here we present an important property that was not 
discussed in [1], namely the particle-particle equal-time 
density-density correlator depending on the distance r: 



Gabir) = 



<vJ:.=i,N.j:j^UNji\rt-r}\-r)> 



(4) 



with ac the Coulomb coupling between charges, Q = E 
or M, and n the density of the corresponding charges. 
It means the ratio of the potential energy (interaction) 



where a, h denotes two species of particles with total num- 
bers Na , Nfj in a volume V and pa , Pb their density. In our 
MD simulation the monopole and anti-monopolc num- 
bers are both set to be 250 inside a sphere with radius 
lOajv/D. a,MD is the length unit in MD and corresponds 
to about 0.2/m after mapping to sQGP system: see [1] 
for more details. The correlator of the monopole-anti- 
monopole pair is particularly instructive because it tells 
how strongly the monopoles are correlated. The result is 
shown in Fig.l, for the 50%-50% electric-magnetic sym- 
metric plasma with F = 2.3. It features a considerable 
nearest-neighbor peak in the monopolc-anti-monopole 
correlator, reaching around 1.25, with hints to small sec- 
ondary correlations: a typical liquid behavior, see various 



examples in e.g. [16]. The same correlator for a weakly 
coupled gas will have a very weak peak, barely above 
1, while a solid will have strong multiple peaks. The 
monopole-monopole correlator shows strong suppression 
at short distance while stays very close to 1 at larger r: 
also typical of Coulomb plasma. The reader shall recall 
that we do not consider here the usual Coulomb plasma 
due to presence of equal number of electric particles in 
our mixture - with which monopoles interacts with a 
strong Lorentz force, leading to very complicated motion 
- so the result is far from trivial. 




FIG. 2. (color online) Monopole-antimonopole (the upper 
two curves) and monopole-monopole (the lower two curves) 
correlators at l.irc(red long dashed) and 3.8rc(blue dot 
dashed): points with error bars are lattice data [3], the dashed 
lines are our fits (see text). 

III. LATTICE DATA AND MONOPOLE 
CORRELATIONS 

Lattice studies of monopoles have a long history which 
we would not even attempt to summarize here (see e.g. 
[17]). We will not discuss here the properties of these 
monopoles - such as masses or their interaction with 
quarks/gluons or "Higgs field" Aq - to be addressed in 
[18] ). We will focus entirely on one aspect of these data, 
related to the correlators (4) and F parameter just dis- 
cussed. Those have been obtained recently in [3], for 
technical reasons for the simplest SU(2) pure gauge the- 
ory. Fig. 2 shows two sets of their data for T = l.lTc and 
T = 3.8Tc for monopole-antimonopole and for monopole- 
monopole correlations. 

Note that the curves are very similar to the ones in 
our MD studies shown in previous section, especially the 
shape and magnitude of the opposite-sign peak. The 
most important feature apparent from these (and other) 
plots is that the correlation gets stronger at the higher 
T, confirming our expectation that the magnetic compo- 
nent of sQGP gets stronger coupled at high T, oppositely 
to the electric component. 



To get a quantitative measure we have fitted them at 
large r (where deviation from 1 is small and linearized 
screening is correct) by the Debye formula 
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FIG. 3. (color online) (a) The magnetic coupling aM 
(on Logio scale) versus T/Tc fitted from the monopole-an- 
timonopole (boxes with solid blue curve) and monopole- 
monopole (triangles with dashed blue curve) correlators. 
Their inverse, the corresponding as from the Dirac condition, 
are shown as stars with solid red curve and diamonds with 
dashed red curve respectively, together with an asymptotic 
freedom (green dotted) curve (see text); (b) Effective mag- 
netic plasma coupling Fa/ versus T/Tc, with the gray band 
schematically showing high-T limit T\.[ (see text). 

In (5) the positive sign in the exponent is for 
monopole-antimonopole {ab = +— ) while the negative 
for monopole-monopole {ab = ++), with um is the mag- 
netic coupling, and the Rd magnetic screening radius. 
See the dashed curves in Fig. 2. Some details about 
the fitting for each data set: we exclude the last few 



points that arc rather flat and marginaUy lower/higher 

than 1 for G_| /G++ which the fitting formulae can 

not reach; the rule to determine how many points from 
large r toward the peaks to be included is to make 
X^/d.o.f as close to 1 as possible; the x^/d-o.f's for 
the eight T's are 1.16,1.96,1.17,0.46,0.82,0.16,1.16,1.23 
for monopole-antimonopole correlators and 3.65,2.87, 
2.16,2.48,2.02,1.06,0.80,1.36 for monopole-monopole cor- 
relators. The values of = obtained arc all in 2—3 
region with slight tendency to decrease toward higher T. 

The resulting values of magnetic coupling from fitting 
both monopole-antimonopole and monopole-monopole 
correlators for all available temperatures are shown in 
Fig. 3 (a): as expected aM is getting weaker at low T 
and stronger at high T. Inverting these values (because 
of the Dirac condition (1)) one gets the respective elec- 
tric couplings, which we compare to the dashed green 
curve corresponding to the pure gauge SU(2) one-loop 
asymptotic freedom expression as ~ (22 / 3) Lnlc*{T/T )] 
with the coefficient C = 54 determined from the last 
point at T = 3.8rc. The one-loop running becomes 
much slower as T decreases, which is normal and shall 
be cured by higher order loops and more importantly by 
non-perturbative corrections such as instanton, see [19]. 

Finally, in Fig. 3(b) we combine the data on the mag- 
netic coupling and the monopolc density from the same 
work to get the dimensionless magnetic plasma parame- 
ter (3) in the studied range T rs l—ATc. We use a formula 
n{T)/T^ = 0.557/[Ln(2.69 * T/Tc)Y from [3] which was 
found to nicely fit the monopole density data. The au 
used is the average of the fitted values from monopole- 
antimonopole and monopole-monopole correlator, with 
half of their difference as the error bar. As one can see, 
magnetic component of QGP never gets to be a weakly 
coupled gas, as Fm > 1 at all T even close to Tc- On the 
other hand, even at the highest T the value of Fm does 
not reach large values > 10 at which liquids are known 
to become glass-like and viscosity starts growing. 

In the high T limit we expect "magnetic scaling" Um ^ 
(q!£;T)^, thus plasma parameter approach a fixed point 

1 /3 

Fm ~ ctMrij^ /T F^ due to Dirac condition. The 
curve in Fig. 3(b) indicate F^ 5. 

IV. SUMMARY 

In this Letter we have shown that gauge theory 
monopoles in a deconfined phase behave as charges in 
a Coulomb plasmas. Furthermore, we show that the 
temperature dependence (running) of the magnetic cou- 
plings in gauge theories is indeed the inverse of the elec- 
tric one, electric-magnetic duality arguments [1]. Good 
agreement was found (in shape and magnitude) of the 
correlators we calculated in MD simulation [1] for novel 
electric-magnetic plasmas and recent lattice results [3]. 
More specifically, we concluded that the magnetic part 
of QGP at T = (1 — 4)Tc has an effective plasma pa- 



rameter in the "good liquid" domain [20] F = 1 — 4.5, 
not spoiling the "perfect liquid" observed at RHIC. We 
predict existence of a fixed point for F at high T. 
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